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The superconductor-insulator transition in the presence of strong compensation of dopants was 
recently realized in La doped YBCO. The compensation of acceptors by donors makes it possi- 
ble to change independently the concentration of holes n and the total concentration of charged 
impurities N. We propose a theory of the superconductor-insulator phase diagram in the (N, n) 
plane. It exhibits interesting new features in the case of strong coupling superconductivity, where 
Cooper pairs are compact, non-overlapping bosons. For compact Cooper pairs the transition oc- 
curs at a significantly higher density than in the case of spatially overlapping pairs. We establish 
the superconductor-insulator phase diagram by studying how the potential of randomly positioned 
charged impurities is screened by holes or by strongly bound Cooper pairs, both in isotropic and 
layered superconductors. In the resulting self-consistent potential the carriers are either delocalized 
or localized, which corresponds to the superconducting or insulating phase, respectively. 



I. INTRODUCTION 

The superconductor-insulator (SI) transition remains 
a challenging and controversial subject after more than 
two decadea 1 ! 2 ' 3 ' 4 ' 5 ' 6 ' 7 ' 8 ' 9 ' 10 ' 11 ' 12 ! 13 . In the low tempera- 
ture limit one can drive the SI transition by changing 
the film thickness, the magnetic field, or the concentra- 
tion of electrons in gated devices. In high T c supercon- 
ductors such as YBCO one can tune the concentration 
of holes by changing the oxygen doping. This leads to 
a SI transition at small hole concentrations of about 6% 
per Cu site in the CuO plane. From this perspective, 
YBCO is essentially a heavily doped semiconductor. It 
is well known that upon decreasing the doping a semicon- 
ductor undergoes the metal-insulator transition when the 
three-dimensional concentration of dopants N crosses the 
threshold Na 3 » 0.02. Here a = h 2 n/me 2 is the effective 
Bohr radius, k is the dielectric constant, m is the effective 
mass and e is the proton charge. In under-doped high T c 
superconductors the conducting phase is a superconduc- 
tor, and one expects a superconductor-insulator transi- 
tion at a similar threshold concentration of dopants. 

In semiconductors one can vary the concentration of 
carriers and impurities independently using compensa- 
tion. For example, in ap-type semiconductor doped with 
Na monovalent donors compensation means addition of 
a concentration Nr> < Na of donors, so that the concen- 
tration of remaining holes n — Na — Nrj becomes much 
smaller than the total concentration of charged impuri- 
ties N — Na + Nd. The metal-insulator transition in the 
(N, n)-plane of a compensated semiconductor was stud- 
ied long ago. It was show n 14 i 15 that in heavily doped 
samples with Na 3 ^> 1 the transition takes place when 
n(N) ~ N/lNa 3 ) 1 / 3 , as was later verified by experiments 
(cf., Fig. 13.3 in Ref.[lJ). 

Recently 1 ^ it was demonstrated that YBCO crystals 
can also be strongly compensated by doping with La. 
Although many of the La 3+ ions substitute for Y 3+ 



and are therefore not electrically active, some La 3+ ions 
substitute for Ba 2+ and hence play the role of mono- 
valent donors compensating oxygen acceptors. It was 
shown that the sample of Yi_ z La z (Bai_ x La x )2Cu30 y 
with x = 0.13 and z = 0.62 is completely compensated 
at y = 6.32, and becomes n-type at y < 6.32. Thus, also 
in high T c superconductors the concentration of impu- 
rities and holes can be varied independently. Resistance 
measurements 1 ^ showed that the SI transition point non- 
trivially depends on both x and y. In strongly compen- 
sated samples it occurs at much larger concentration of 
holes than in standard uncompensated samples. How- 
ever, the full phase diagram of the zero-temperature SI 
transition in the plane (N, n) has not been established 
yet experimentally. In this paper we predict it theoreti- 
cally. 



A. Global phase diagram 

Let us start by discussing the gross features of the 
phase diagram which are expected, e.g., in compensated 
high T c materials such as La doped YBCO, see Fig. [T] 
In the uncompensated material with n — N, we expect 
a transition from the insulator to a superconductor at a 
critical doping n u (on the underdoped side) as discussed 
above. The pairing mechanism is believed to be at least 
in part due to spin fluctuations which become signifi- 
cantly weaker upon exceeding an optimal doping level. 
Finally, superconductivity is essentially destroyed on the 
overdoped side (n > n ), or at least T c is strongly sup- 
pressed. Upon adding the temperature axis to the phase 
diagram this leads to the well-known superconducting 
dome in high temperature superconductors. As disorder 
is increased by compensation (increasing N/n), the dop- 
ing concentration n u (N) where delocalized states first 
appear, increases as well. On the other hand, we expect 
that the upper critical density n (N) decreases because 
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usually disorder diminishes the effectiveness of the super- 
conductive attraction, while it enhances the competing 
Coulomb repulsion. We thus propose that at some com- 
pensation N/n — N* Jn where n u (N*) — n (N*) = n* 
there may exist a tricritical point beyond which a direct 
transition from a localized insulator to a metal without 
intermediate superconducting state takes place. Note 
that the effect of compensation is similar to that of a 
strong magnetic field: both suppress superconductivity. 

In this paper we are not concerned with the transi- 
tion to a metal at high doping, nor with the vicinity 
of the tentative tricritical point T = {n u (N*), N*} in 
Fig. [TJ where metal, insulator and superconductor meet. 
Instead we analyze the dependence of n u on the degree of 
compensation. In the case of strong coupling supercon- 
ductivity the latter exhibits interesting new features in 
the regime of low densities, reflecting the crossover from 
a BEC (Bose-Einstein condensate) to a BCS supercon- 
ductor in the interacting gas of preformed Cooper pairs. 



B. BEC-to-BCS crossover in the SI transition 

A first attempt to predict the low density part of 
the SI phase diagram^ was based on the toy model of 
an isotropic compensated p-type semiconductor with a 
strong (unspecified) pair-forming mechanism. The size 
of hole pairs £ was taken as a free parameter as deter- 
mined by a strong coupling mechanism. For the major 
part of this paper we will adopt this approach as well. For 




FIG. 1: Global phase diagram of compensated high T c su- 
perconductors in the plane (n,N). Here M, S and I stand 
for metal, superconductor and insulator, respectively. The 
line n u (N) separates the insulator form a superconductor, 
which eventually turns into a metal upon overdoping beyond 
n (N). We conjecture the existence of a tricritical point 
T = {n u (N*),N*}, where n u (N*) = n a (N*) beyond which 
the insulator turns directly into a metal. A large part of the 
paper is focused on the lower left corner of the line n u {N) 
which exhibits features of the BEC-to-BCS crossover of the 
interacting gas of Cooper pairs in strong coupling supercon- 
ductors. 



simplicity we assume £ to be independent of the density 
of carriers and the disorder, at least in the dilute BEC 
part of the phase diagram. However, it would not be dif- 
ficult to account for such a dependence (in certain strong 
coupling models for preformed pairs one expects such a 
dependence even in the dilute regime^) . In the following 
we thus concentrate on the two independent variables n 
and N, taking £ as a fixed parameter. 

In Ref. [I?] two limiting cases of the SI transition were 
identified: In the limit of large pairs which overlap sig- 
nificantly in space, n£ 3 3> 1, one obtains the stan- 
dard Bardeen-Cooper-Schrieffer (BCS) instability of the 
fermion system. If disorder is weak the electrons are de- 
localized and form a dirty BCS superconductor. This 
happens essentially at the same critical density as the 
metal-insulator transition in a semiconducto r 14 ! 15 with- 
out superconductivity, 

iV 

Note that there is no dependence on £ in {I), because 
electrons are only weakly bound and, therefore, screen 
the random potential of charged impurities like free ones. 
Here and in all formulae below we omit numerical coef- 
ficients and adopt the scaling approach. The scaling is 
controlled by the large dimensionless parameter Na 3 3> 1 
and the dimensionless ratio a/£. 

The opposite limit of very small and strongly bound 
pairs is more unusual. Upon decreasing the concentration 
of holes n the SI transition occurs due to the localization 
of hole pairs (composite bosons) in a random potential 1 . 
At small external disorder the bosons undergo a Bose- 
Einstein condensation (BEC) while at large disorder the 
condensate is fragmented and turns into a Bose insulator 
(also referred to as Bose glass). This limit is reached 
when the pairs are dilute, n£ 3 <C 1, and can be considered 
as a gas of point-like charged bosons. A similar picture 
applies to the case of neutral boson o 18 i 19 . 

As we will rederive below, for point like bosons the 
border between the superconducting BEC phase and the 
Bose insulator occurs at the hole density 

n = n >w=(Nwrs> (2) 

where typical screened Coulomb wells loose their bound 
quantum levels. Of course, the length £ is again ir- 
relevant, because the pairs are considered as point- like 
bosons. Notice that for a heavily doped system with 
Na 3, 3> 1, we have na(iV) 3> n\(N). In other words, in 
a given disorder, a system of small pairs delocalizes at a 
much higher density n than a system of weakly bound 
electrons with larger pair size £. This reflects the fact 
that, at equal density n, bosons have less kinetic energy, 
and thus one needs more of them to induce their collec- 
tive derealization. 

The crossover between the above two limiting cases 
is quite subtle. In Ref. [l?] it was incorrectly assumed 



3 



na J 



n A a-- 



1 



n > . 



n, 



A. 



B 



n 2 



I 



n^ 3 = 1 



1 N B a3 



N A a 3 



Na 3 



FIG. 2: Si-phase diagram in an isotropic 3D system, on a 
log- log plot (S stands for superconductor, I for insulator). 
The dilute boson part of the curve (C-B) is described by nz, 
Eq. the interacting boson part (B-A) by ri2, Eq. (|19[) . and 
the standard strong coupling BCS transition beyond A by ni, 
Eq. dXJ) . The BEC part of the transition line (C-A) only exists 
if the pairs are very small, £ < a. In the case of larger pairs 
(£ > a), and for weak coupling in general, the ni line extends 
all the way to the point C. The shaded region corresponds to 
n > N which is unphysical. 



that the limits described by rii(N) and ns(iV) require 
only the inequalities n£ 3 > 1 or n^ 3 < 1 to be valid, 
and hold all the way up to the BEC-BCS crossover line 
rt£ 3 = 1. However, this argument neglected the repulsion 
between bosons as they become denser, and thus it lead 
to the incorrect conclusion that the BEC-BCS crossover 
line forms a substantial part of the SI transition line. 
Below we reconsider this crossover in detail. 

This paper contains two main new results. First, in 
Sec. [IT] we show that while rii(N) is valid all the way up 
to the BEC-BCS crossover line < 3 = 1, the BEC part 
of the SI transition is to a large extent dominated by an 
intermediate segment n = n-z (N) of the transition line at 
which the chemical potential of repulsive compact bosons 
becomes of the order of the amplitude of the screened 
Coulomb potential. This segment interpolates between 
the above discussed limits ni(N) and ri3(7V), see Fig. [2] 
We will see that the BEC regime, ri2(N) and ri3(N), 
occurs only when the pairs are smaller than the effective 
Bohr radius, £ < a. We show below that in this case 
transport in the insulating phase is due to the hopping 
of hole pairs. In App. [X] we will discuss bipolarons as an 
example of strong coupling superconductivity which can 
give rise to such small pairs. 

Second, in Sec. IIIII we apply similar ideas to a 
generic strongly anisotropic superconductors with a lay- 
ered structure, such as formed by the CuO or FeAs planes 
in high T c superconductors. We arrive at a qualitatively 
similar phase diagram in the plane (N, n) for this case 
as well, see Fig. @] The details of the phase diagram are 
found to depend on the ratio between Bohr radius and 
interlayer distance. 



II. SI PHASE DIAGRAM OF AN ISOTROPIC 
SUPERCONDUCTOR 

Let us now recall the derivation of the limiting critical 
concentrations ni(N) and n,3(N). We consider the case 
of heavily doped materials, Na 3 1, which provides 
a large parameter that makes the scaling analysis well 
controlled. 



BCS segment of the superconductor-insulator 
transition line 



We start from the BCS side at high density (large 
pairs). Let us divide the sample into cubes of linear size 
R. Due to spatial fluctuations of the concentrations of 
donors and acceptors each cube contains a random impu- 
rity charge of arbitrary sign and with an absolute value 
of the order of e(N i? 3 ) 1 / 2 . At the scale R such randomly 
fluctuating charges create a random potential energy re- 
lief of amplitude 



zV{R) 



e 2 {NR 3 ) 1 / 2 _ e 2 {NRf' 2 



kR 



(3) 



This energy diverges at large R, so that screening even 
by a small concentration of holes n is crucial. To discuss 
this screening we estimate the characteristic fluctuating 
density 5N(R) of impurity charges at the scale scale R: 



SN(R) 



(NR 



3U/2 



R 3 



N 
R 3 



1/2 



(4) 



The concentration n of carriers can be redistributed be- 
tween wells and hills of the random potential. This redis- 
tribution screens all the scales R for which SN(R) < n 
or, in other words, for R> R s , where 



1/3 



(5) 



is the nonlinear screening radius^±£. All scales R < R s 
remain unscreened, because even when all electrons are 
transferred from all the hills of the potential energy to all 
its wells they are not able to level off the charge density of 
such fluctuations. Since V(R) cx i? 1 / 2 among remaining 
scales the most important contribution to the random 
potential is given by R = R s . Thus, the amplitude of the 
nonlinearly screened random potential energy is 



eV{R s ) 



e 2 N 2 ' 3 

K 77, 1 / 3 



(G) 



So far we have dealt only with the electrostatic energy 
of holes and neglected their kinetic energy. At T = 
all kinetic energy is of quantum origin, and we should 
find the conditions under which it is small enough so 
that the above described picture of localized electrons is 
valid. Clearly the potential energy Eq. ^ is able to lo- 
calize electrons with concentration n if it is larger than 
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the Fermi energy of holes in its wells ep(n) = h 2 n 2 ^ 3 /2m 
(m being the effective mass). In the opposite case 
ep(n) 3> eV(R s ) the Fermi sea covers the typical maxima 
of the potential energy relief and the semiconductor be- 
haves like a good conductor, see Fig. [3] Equating eV{R s ) 
and e F (n) we obtain the critical concentration ni(N) for 
the SI transition, as given by Eq. (fT]) 14 i 15 . Note that the 
non-linear screening theory requires nR 3 — N/n 3> 1 
which is always fulfilled in strongly compensated materi- 
als. 

In the delocalized phase electron screening becomes 
linear, the screening radius being given by the standard 
Thomas-Fermi expression 



1 d[i\ 
e 2 dn J 



1/2 



(na 3 ) 1 / 6 ' 



(7) 



and the amplitude of the screened potential relief equals 
eV{r s ) — e 2 (Nr s ) 1 / 2 / n. As expected, at the transition 
n = ni, these two quantities match the corresponding ex- 
pressions Eq. ([5]) and Eq. ^ pertaining to the insulating 
side. 



BEC segments of the superconductor-insulator 
transition line 



In the above discussion the notion of strong pairing 
attraction and preformed pairs was irrelevant. However, 
as we follow the transition line Eq. (fTJ) to lower densities, 
we may finally reach the crossover to the BEC regime, 
which takes place when strongly bound hole pairs be- 
come dilute, i.e., when n^ 3 = (Na 3 ) 2 ^ 3 (^/a) 3 = 1. This 
corresponds to the point A in Fig. [2] and the densities 



(a/0 3 , 

(«/0 9/2 . 



(8) 
(9) 



Under the assumption of heavy doping, Na 3 1, 
the crossover to the BEC regime can only happen when 
the pair size is much smaller than the Bohr radius, 
£ < a(Na 3 )~ 2 / 9 < a. For the sequel we will assume that 
the pairs are very small S, < a. Since we will be using the 
concept of strongly bound pairs a lot, we briefly recall the 
essential elements of strong coupling superconductivity. 



1. Strong coupling superconductivity 

The physics of a fermion gas subject to attractive in- 
teractions (but in the absence of disorder) has been stud- 
ied in detail in Refs. [20II21II22I and is now a very active 
field of studies in the context of cold atoms 23 . The au- 
thors of Ref. considered electrons with a mutually 
attractive potential of size Vk,k' ~ V for k < ko, and 
rapidly decaying for larger k. If the interaction poten- 
tial between two holes is too weak to produce a bound 
state (V < V c ~ 1/mko), the fermions are essentially 



unbound, and only an exponentially narrow range of en- 
ergies around the Fermi level participates in pairing, the 
gap being of the order of 



A 
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-^exp{-l/v(E F )V kj 
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(10) 



where v(Ep) is the density of states at the Fermi level. 
On the other hand, if the mutual interaction is strong, 
bound states of two single carriers exist, and at low den- 
sity the fermions organize into preformed pairs with a 
typical size £ » [{V/V c — l)^] -1 and a pairing energy 
Spair ~ h 2 /m{; 2 . As long as the pairs are dilute, n£ 3 <C 1, 
the chemical potential for the addition of pairs, fi, and 
the gap function A are much smaller than the pairing 
energy, 



A 



ft 2 

£pair « 3 ) n£, 

m 

£pair K 3 ) 1/2 - 



(11) 
(12) 



However, when the pairs become dense, the pair chem- 
ical potential is dominated by the Fermi energy of its 
constituting fermions, 



/1 ~ Ep = 



h 2 n 2 ' 3 



rn 



(13) 



At the same time the gap function A crosses over to its 
strong coupling BCS form, i.e., Eq. (|10p with an exponent 
of order^i O(l), 



E F . 



(14) 



In this dense regime the pairing energy is dominated by 
the gap function -E pa i r ~ A > h 2 /m£ 2 . 



2. Very dilute bosons 

In order to derive the critical concentration 713 (N) of 
Cooper pairs (charge 2e bosons) at the SI transition, we 
notice that the above calculation of the nonlinear screen- 
ing radius R s ([5]) and the random potential energy cre- 
ated by screened charged impurities © remains unal- 
tered in the scaling sense. 

The difference between the gas of composite bosons 
and that of weakly bound fermions lies in their quantum 
kinetic energy^. Due to the weak effect of Pauli's prin- 
ciple on strongly bound Cooper pairs, a large number of 
them can occupy a given localized level of a potential 
well, keeping the quantum kinetic energy low. There- 
fore, the condition of de-localization of Cooper pairs is 
much more stringent than the condition eV(R s ) < ep (n) 
which applies to fermions. A sufficient condition for the 
derealization of a compact Cooper pair is that a typical 
well of the random potential does not contain any local- 
ized level, or eV{R s ) < h 2 /mR 2 , where m is the effective 
mass of pairs, which we assume to be of the same order 
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as that of electrons. This condition is also necessary if 
mutual repulsions can be neglected, as we will discuss 
below. Solving the equation 



eV(R s 



mR 2 



(15) 



for n and using Eqs. ([5]) and ([6]) we find the critical 
concentration of the SI transition given in Eq. $Z§. This 
derivation clearly demonstrates why n${N) ^S> ni(N). 
According to Eq. ^ the potential energy amplitude 
eV(R s ) decreases with increasing n. To achieve der- 
ealization it has to be pushed below the quantum kinetic 
energy of the clean system. This requires larger n in the 
boson case and thus leads to nz(N) ^> ni(N). 



3. Moderately dilute, interacting bosons 

So far, following Ref.ll7L we have taken into account all 
the Coulomb interactions. However, we have neglected 
the short range repulsive interaction between composite 
bosons. Such a repulsion is related to the Fermi nature 
of individual holes, which becomes important if two pairs 
of holes overlap within their length £. This short range 
interaction can be described by the well-known expres- 
sion for the chemical potential fi(n) of a non- ideal gas of 
bosons of concentration n with a scattering length £: 



ji = (h 2 /m)n£, 



(16) 



which is also confirmed by the result (fTTj) for dilute sys- 
tems of strong coupling superconductivity. This chemical 
potential reflects the extra quantum kinetic energy due 
to the mutual repulsion of the pairs. Note that it matches 
the Fermi energy E F = H 2 n 2 / 3 /m when the BEC-BCS 
crossover n£ 3 = 1 is reached. 

The derealization criterion (p~5|) discussed above re- 
mains relevant as long as the density is low enough such 
that fi is smaller than the typical localization energy 
Ti 2 /mR 2 . However, at an impurity density N — Nb, the 
chemical potential of the critical insulator (n = 713 (N)) 
becomes of the order of the typical amplitude of the ran- 
dom Coulomb potential Eq. (|6|). This marks the point 
B in Fig. [51 beyond which the derealization is driven by 
the mutual repulsion between bosons. The crossover in 
the transition line occurs at the densities 
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B 



(a/0 2 , 
(a/0 



5/2 



(17) 
(18) 



On the low density side the ri3(N) line ends at point C 
which corresponds to the uncompensated limit 71.3a 3 = 
Na 3 — 1. At higher densities, N > Nb we need to 
compare the quantum kinetic energy Eq. (fT6|) to the am- 
plitude of potential fluctuations, Eq. ©, similarly as in 
the BCS regime. This leads to the new segment of the 
transition line 



n 2 (N) 



N 1 ' 2 

K) 3/4: 



(19) 



which interpolates between points A and B in Fig. [2j 

We can confirm this result by calculating the linear 
screening radius r s in the delocalized Bose gas. Us- 
ing Eq. (|16p to compute the compressibility we find the 
Thomas-Fermi screening radius from Eq. j7} as r s = 
(a^) 1 / 2 . One easily verifies that this linear screening ra- 
dius matches the non-linear screening radius ([5]) at the 
transition line (TH))) . Similarly, one can check that along 
the line n\ the linear screening radius of the conducting 
side matches the non-linear screening radius R s on the 
insulating side. The linear screening in the very dilute 
superfluid above the line 713 is found^ to be 



n 1 



(20) 



which again matches R s at the transition line 713. 

The fact that the SI transition line undergoes a kink 
at the BCS-BEC crossover is very similar to the case of 
the superfluid-insulator phase transition in a neutral gas 
of attractive fermionsi^. However, at lower density 712 
lies below the BCS-BEC crossover line, contrary to what 
was claimed in Ref. 

The results obtained so far in this section can be sum- 
marized in the following concise manner: The chemical 
potential of a gas of composite bosons of size £ , localized 
into a region of linear size R is given by 



/i(n, R) 



iR 2 ' 



-nmin(£,n 1/3 ) 



(21) 



The first term of the righthand side refers to the ground 
state energy in a well of size R. At higher density /j,(n, R) 
is dominated by the second one, describing the interac- 
tion energy §Tty of repulsive bosons (for n£ 3 < 1), and the 
Fermi energy Ep ~ n 2 / 3 of the BCS regime (for n£ 3 > 1), 
respectively. 

The SI transition occurs when the chemical potential 
dominates over the amplitude of the screened impurity 
potential, i.e., when 



fi{n,R s )~eV(R s ). 



This can be reformulated as 



n 3 (N)J ' m(A0 



i[i,« 3 K 



(22) 



1, (23) 



which defines the transition line in the whole (n, ./V)- 
plane, as plotted in Fig. [51 



C. The nature of the insulating regime 

1. Droplets in the insulator- 
It is important to understand the insulating phase in 
some more detail. Deep in the insulator, the charge den- 
sity n is by no means homogeneously distributed. Instead 
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the holes or Cooper pairs fill deep wells, where they form 
puddles of high density, while the rest of the space is 
completely void of carriers. 

To determine the chemical potential and the size of 
puddles we can argue as follows: when n <C n(N) only 
the deep wells of the landscape are populated with car- 
riers (see Fig. [3J . 

Suppose the carriers fill a well of linear size R < R s . Its 
typical depth eV(R) is given by Eq. (J3J) , and it contains 
an excess impurity charge of order eQ(R) = e{NR 3 ) 1 / 2 . 
Upon filling the well with carriers, their chemical poten- 
tial raises continuously with respect to the bottom of the 
well. Assume that when the chemical potential reaches 
eV(R) we have filled in Q q (R) particles. If R is small, 
Qq(R) < Q{R), that is, the exclusion principle or the re- 
pulsion between bosons limits the number of carriers we 
can fill into the well. In large wells, we can at most fill in 
Q(R) particles before turning the well into a hump. How- 
ever, those larger wells will not be filled homogeneously 
with particles. Rather, they split into smaller droplets 
for which Q q (R) w Q(R). The latter relation define: 
typical droplet size R q : 



In the repulsive boson regime, (A£ 3 ) < n£ 3 < 1 (or 



,15 



A* Tl 



Q[Rq 

R q 



= eV{R q ). 



(24) 



Here we have to use the expressions for the chemical po- 
tential given above in Eq. (j2~l]) . 

Analyzing the three regimes of the phase diagram 
(BCS side, repulsive bosons and very dilute bosons) we 

find the following results. When n£ 3 > 1 (A > on 
the transition line), fj, is given by Ep(n), and we find the 
droplet sizei 5 . 



R q ,i — 



(a 3 7V)V9- 



(25) 





r\ / 






R 


q 


eV(R s ) 




FIG. 3: Top: Droplets in the insulator. The carriers assemble 
in non-percolating fragments of size R q < R s . Bottom: In 
the conductor, the wells are not deep enough to localize the 
carriers. The latter delocalize due to their quantum kinetic 
energy, n > V(r s ) (fj. = Ef in the BCS limit of dense pairs). 



M (3d) 

stead 



< A < A 



(3d) 



on the transition line), we find in- 



R 
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(26) 



Eventually, in the lowest density regime, n£ 3 < (A£ 3 ) 2 
where the bosons do not significantly interact, R q is sim- 
ply the typical localization radius in the disorder poten- 
tial. It is obtained from h 2 /mR 2 ~ eV(R) as 



R 



9,3 



(a 3 AO 1 / 5 ' 



(27) 



In all three cases, the insulator consists mostly of pud- 
dles of size R q which are well separated and do not per- 
colate. One can verify that the SI transition occurs when 
the droplets grow to the size of the non-linear screening 
radius, R q = R s = (A/n 2 ) 1 / 3 . Indeed, at this point 
droplets of size R s start to percolate, which induces the 
derealization transition. 

Note the remarkable fact that in all insulating regimes 
the density of carriers in the above droplets is the same 
as the critical density of the corresponding segment of 
the SI transition, n q = Q{R q )/R% = n(N). 



2. Level spacing in droplets 

Before we turn to the role which droplets play in the 
transport properties of the insulator phase, we have to 
discuss the level spacing in a typical droplet. The typ- 
ical cost to add another carrier into a droplet is S = 
R~ 3 d/j,/dn(n = n q ). This is essentially the level spac- 
ing of the considered droplet. Interestingly this quantity 
turns out to be equal to the charging energy e 2 /R q . This 
holds both in the dense BCS-like part and in the inter- 
acting boson regime of the phase diagram. In the very 
dilute boson regime, the quantity of interest is not the 
level spacing S, but the typical kinetic energy scale of 
a localized wavefunction, which again turns out to be 
equal to the charging energy e 2 /R q . Thus, for the above 
(spontaneously originating droplets) we do not have to 
distinguish between one-electron level spacing and charg- 
ing energy. 

This unique energy scale is important in determining 
whether all carriers are paired or whether it is energet- 
ically favorable to break up a pair and redistribute the 
constituting holes onto two different droplets. The cost 
of such a break up is the pairing energy, while the maxi- 
mal energy gain is of order S ~ e 2 / ' R q . Thus the criterion 
for having all carriers paired up in the ground state is 



2 /R q < E pa 



(28) 



In the BEC regime the pairing energy is given by 
-Epair ~ h 2 /m£, 2 . In the BCS regime it is even bigger, 
if the BCS coupling remains strong, see Eq. (|14p . Thus, 
if £ < a, it is never favorable to break Cooper pairs, i.e. 
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the insulating state is always a Bose glass. Moreover the 
droplets are actually superconducting at low enough tem- 
peratures. In tunneling experiments they should show 
a hard gap with coherence peaks on its shoulders, de- 
spite the absence of global phase coherence among the 
droplets. 



D. Variable range hopping transport in the 
insulator 



The above implies that for systems with small pairs, 
£ < a, the low temperature transport is Efros-Shklovskii 
variable range hopping of Cooper pairs between droplets. 
This yields a conductivity 



<t(T) = ctq exp 



1/2 



with a characteristic temperature 

(2e) 2 



Tes = 2. 



k£ 2 



(29) 



(30) 



where l 2 is the effective localization length of a Cooper 
pair. This prediction agrees qualitatively with experi- 
mental datai&2&. 

However, there is an exception to the above assertion 
that pairs prevail in the insulator if £ < a. Namely, if the 
dimensionless BCS coupling A decreases with increasing 
density in the BCS regime, the pairing energy -E pa i r can 
become exponentially suppressed at high densities. (This 
presumably happens on the overdoped side of high T c 
superconductors.) 

When -Ep a i r ~ A falls below the level spacing in typical 
droplets, 5, it becomes favorable to break up pairs and 
redistribute the carriers on different droplets. One can 
verify that at the same time the parity gap (the extra 
cost for having an odd number of particles per droplet) 
becomes smaller than the level spacing 2 ^. In this situa- 
tion the ground state of the system is a Coulomb glass of 
unpaired fermions. Consequently, the low temperature 
transport is again of the form (|2"5|) , but with a character- 
istic temperature 



T ES = 2. 



(31) 



which is roughly 8 times smaller than (|30|) . because the 
localization length of a hole is about twice as big as that 
of a pair, t\ ~ 2i 2 - When pairs are not very strongly 
bound (£ > a), as well as in the case of weak coupling, 
only the BCS segment of the superconductor-insulator 
border survives. In the weak coupling case, it can easily 
happen that £/ pa ir < 5 in the droplets of the insulator. In 
this case they contain odd or even numbers of holes, and 
the low temperature variable range hopping is dominated 
by unpaired holes. 



So far we have not specified any particular strong cou- 
pling mechanism which leads to the preformed pairs on 
the insulating side of the SI transition. In App.[A]we dis- 
cuss an explicit example of strong coupling superconduc- 
tivity, which allows one to formulate a direct microscopic 
criterion for the condition £ < a. As discussed above, 
the latter is necessary to observe the BEC part of the SI 
transition in a heavily doped system. 



E. Coulomb correlation energy 

As we saw above the long range Coulomb interactions 
considered within mean field approximation play a ma- 
jor role in our theory. However, the correlation energy 
produced by the Coulomb interaction between nearest 
neighbors has been neglected so far. We should thus 
make sure that corrections to the chemical potential due 
to Coulomb correlations are subdominant with respect to 
the leading term given in Eq. dHJ). According to Foldy 2 ^ 
the energy per particle in a disorder-free, Coulomb inter- 
acting Bose system in 3d is 



= f_ p7 3/4 = !L n 2/3 p 5/8 = ^ n l/4 a -5/ 4j 



(32) 



where p s = (a 3 n) -1 / 3 . One can verify that this quantity 
is indeed smaller than /x(n, R s ), along the whole SI phase 
transition line, consisting of the segments ni(N), n 2 (N) 
and n 3 (N). 



III. SI PHASE DIAGRAM OF A LAYERED 
SUPERCONDUCTOR 

A. General theory 

1. Non-linear screening in a layered system 

In this section, we extend our previous arguments to 
the case of anisotropic layered superconductors. We as- 
sume that with respect to the motion along the c-axis 
(z-axis) all holes reside in the lowest spatial quantiza- 
tion mode of the narrow quantum wells defining the 
conducting a6-planes (x, y-plane) perpendicular to the 
c-axis. These parallel wells are located at a distance d 
from each other, each well containing holes with the two- 
dimensional concentration nd. Impurities of both signs 
are randomly distributed between these narrow quantum 
wells. We assume again that a strong attraction between 
the holes of a given well leads to preformed pairs (com- 
posite bosons) with a size £ in the plane of the well. To 
simplify things we will first assume an isotropic dielec- 
tric constant k. Modifications due to anisotropy will be 
discussed in the next subsection. 

Let us define again the effective Bohr radius in the 
a6-plane, a = h 2 K/me 2 . For the insulating phase we 
need to understand the non-linear screening in a system 
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containing impurities in the bulk and screening carriers 
confined to planes. There are two limits of this screening 
problem. When the non-linear screening radius is bigger 
than the distance between layers, R s > d, we can use our 
results for the isotropic 3d case, cf. Eqs. (|5I6|) . On the 
other hand if d > R B , the potential fluctuations within 
each plane are screened independently. For these two 
cases we have obtained the phase diagrams of Fig. 2] 
The specific expressions for the various lines are derived 
below. 

The non-linear screening radius for the case d > R s 
has been derived in Refs. l3dl)3lL Let us cover a conduct- 
ing plane by densely packed cubes of linear size R < d. 
Fluctuations of the charge among these cubes are of the 
order of e(NR 3 ) 1 ^ 2 . The random potential they create in 
the planes can be screened by redistributing the charge 
of two-dimensional holes, endR 2 between potential hills 
and wells of linear size R, if the latter is large enough. 
We find the nonlinear screening radius R s by equating 



e(NR 



3U/2 



n^ 2 d = 1 




l^a J N B a J N A a 



Na 3 



n^ 2 d = 1 




FIG. 4: The superconductor-insulator transition in the plane 
(n, TV) for the model of a layered high T c superconductor (on 
a log- log scale). S stands for superconductor, I for insulator. 
The shaded region n > TV is unphysical. (a) Narrowly spaced 
layers, d < a. At lowest densities the SI transition is a Mott 
transition in the layers. The very dilute boson part of the 
curve (M-B) is described by 715, Eq. (|37[1 . while the interacting 
boson and BCS parts (B-A and beyond) are described by ru, 
Eq. (|36p . (b) Widely spaced layers, d > a. The dilute boson 
regime does not exist. In both cases, the point A where a weak 
BEC-BCS crossover takes place, exists only if (ad) 1 / 2 > £ or 
a > f, respectively. 



endR 2 , which yields 
TV 



(nd) 2 ■ 



(33) 



Only scales of the random potential with R < R s sur- 
vive the screening. Thus, the amplitude of the remaining 
random potential is 



K nd 



(34) 



At small enough hole concentration n the screening ra- 
dius is always bigger than d. However, upon approaching 
the SI transition, either of the above screening scenarios 
may apply. As we will see, the first case applies to small 
separations of layers, d < a, while the screening of inde- 
pendent layers governs close to the SI transition if d > a. 



2. Narrowly spaced layers, d < a 

The difference between layered materials and isotropic 
ones is due to the confinement of the carriers to the layers. 
The quantum kinetic energy, or the chemical potential of 
pairs confined to a region of linear size X in the plane, 
evaluates to 



fi(n,X) 



nd 



1X 2 ' to max [1, log(l/£ 2 nd)] 



(35) 



Note that in 2d there is only a logarithmic difference 
between the Fermi energy Ep in the fermion fluid and the 
interaction energy per particle in the Bose gas, both scal- 
ing essentially as (H 2 /m)nd. The logarithmic dependence 
of on the pair size £ in the interacting boson regime is 
well-known^ 2 .. Note that the logarithm is replaced by 
unity at the BEC-BCS crossover point nd£ 2 = 1. The 
first term in (|35[) is the kinetic energy of a single boson 
or fermion in a well of size X, the relevant size in the 
context of non-linear screening being X = R s . 

We will show below that the logarithmic effects due to 
the finite size of pairs £ is only relevant for the transi- 
tion when the pairs are small, £ < (ad) 1 / 2 . Let us thus 
first discuss the opposite case £ > (ad) 1 / 2 . The chemical 
potential for pairs, Eq. (|35p . then scales in the same way 
as that for fermions without any superconducting corre- 
lations. Thus the results below describe equally well the 
metal-insulator transition of unpaired fermions confined 
to planes. 

Delocalization and hence the insulator-conductor tran- 
sition, takes place roughly when u(n, R s ) — eV s . In the 
high density regime where the Fermi energy dominates, 
the transition occurs at 



n 4 



TV 1 / 2 
(ad) 3 /4 ' 



(36) 



The difference between this result and Eq. (fTJ) is due 
to the confinement of holes to the 06-planes. At low den- 
sities, the random potential eV s competes against the ki- 
netic energy due to the confinement of carriers to regions 
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of size R s in the plane. ^From this we find the critical 
density 



N 



n 5 



(TVa 3 ) 1 ^' 



(37) 



which is the same as in the isotropic case, Eq. ©. 

The line (J5J) continues down to densities Um = l/cs 2 d 
which is the minimal density required to drive a Mott 
transition in the layers. Note that for d < a this minimal 
density is higher than in the isotropic case, since here the 
carriers are confined to narrowly spaced planes. 

The crossover between ri4(N) and ns(AT) occurs at the 
densities 



(2d) 3 

n y B 'a 



N {2d) a 3 



■V 2 



(38) 



The transition lines 774 and 775 can also be derived by 
approaching from the conducting side, in close analogy 
to the isotropic case. 

It is justified to deal with fermions and thus to ignore 
the logarithm in Eq. (|35|) if we are still on the BCS side 
of the BCS-BEC crossover, i.e., if ^ 2 nd > 1 down to 
ri = rig . This condition is equivalent to £ > (ad) 1 / 2 as 
we anticipated above. In this case, the crossover to the 
dilute boson (BEC) regime occurs along the line 715 only, 
without affecting the shape of the SI transition line. 

Let us now discuss effects which arise if fermions bind 
into small pairs of size £ < (ad) 1 / 2 . As in the isotropic 
case, there is a BCS-BEC crossover at the point A in the 
phase diagram, where n^ 1 = 1/ £ 2 d > n^' and 



(2d) 3 

n A a 



ed 1 



(3) 



1/2 



(39) 



However, here the difference between the fermion regime 
and the interacting boson regime results only in a loga- 
rithmic factor correcting the line 714 to 



714 = 



(ad) 3 / 4 S yedN 1 / 2 J 



(40) 



3. Widely spaced layers, d > a 

If the spacing between layers is larger than the Bohr ra- 
dius we need to compare the potential fluctuations (|3"4"|) 
to the chemical potential ([35]) . As above, we will find 
that when pairs are large, £ > a, they do not affect the 
phase transition line, which then becomes equivalent to 
the metal-insulator transition of an unpaired fermion sys- 
tem. In the high density regime, equating Ep to eV s from 
Eq. (j3~4"|) , we find the transition line 



774 



N 1 ' 2 
aM 2 d 



(41) 



The non-linear screening radius remains constant R s = a 
along the SI transition line. 



It turns out that, contrary to the case d < a discussed 
above, the line (1411) describes the SI transition down to 



dopant densities where Nm& = 1 and 77 = 77m = X/a d, 



while the first term in (|35[) never becomes relevant. Once 
the dopants are dilute, Na 3 < 1, one leaves the regime 
of heavy doping. The potential for individual carriers is 
then dominated by the closest impurity charge. Under 
these conditions the derealization takes place as a stan- 
dard Mott transition in the planes. It occurs when the 
nearest neighbor distance in the planes is of order a, i.e., 
when 



1 



(42) 



Again, it is justified to deal with fermions and to neglect 
logarithmic factors if £ 2 nd > 1 holds down to 77 = um- 
This is equivalent to the condition of large pairs, £ > a. 

However, in the case where fermions are strongly 
bound into pairs of size £ < a, there are logarithmic 
corrections to the phase boundaries. One finds that the 
line (141[) turns slightly upward, 



ATl/2 



774 



log 



r l/2 



7 1 / 2 d 



(43) 



C ^ (2d) 

lor 77m < n < n A 
shown in Fig. UJd. 



a i/2 d * y^2 dN i/2^ 

l/£ 2 d. The full phase diagram is 



B. Anisotropic dielectric constant 

Here we refine the above analysis and take into account 
the anisotropy of the dielectric constant in a layered sys- 
tem. We use k z for k zz and k x for k xx — n yy . We also 
define the average dielectric constant as k — (k^Kz) 1 . 

In order to derive the SI transition line in the presence 
of an anisotropic dielectric constant we first switch to the 
new coordinate frame (2/ , y' , z'), where x' — xj ' n]J 2 , y' — 

yl^J 2 and z' = x/kI^ 2 . In this framed the Coulomb 
interaction of a charged impurity with a hole becomes 
isotropic e 2 /^ 3 / 2 /. At the same time, the concentrations 
N and 77 are transformed, too: N' = k 3 / 2 N, n' = K 3 / 2 n. 

Let us first treat the case of narrowly layered sys- 
tems. In the new frame we have similarly to Eqs. (|5l6p 

R' s = AT'V3/ n /2/3 = ^1/3/^1/2^2/3 and eV (R! a ) = 
e 2jV' 2 /3/ K 3/2 n /l/3 = g 2 iv 2/3/ Kn l/3^ Thuf ^ returnm g to 

the laboratory system we arrive back at Eq. (J5J) for the 
amplitude of the screened potential eV s , and to the same 
SI transition line (|36| with redefined k = (k 2 /-^) 1 / 3 . 
However, note that the notion of the nonlinear screening 
radius R s becomes anisotropic. Characteristic potential 
wells have a scale X s = k x /2 R' s = (N 1 / 3 /n 2 ' 3 )a l l G in 
the (x,y) plane, where a — k x /k z . On the other hand, 



the scale perpendicular to the planes is Z s 



1/2 



(jyi/3/ n 2/3) a -i/3 = x s a- x l 2 < X s . This anisotropy 
modifies the critical concentration 775 (i.e., the very di- 
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lute boson limit) to 



ri5 



(Na 3 ) 1 / 5 ' 



(44) 



It also affects the criterion on the smallness of pairs that 
is required for a BEC regime with logarithmic factors 
to exist. The criterion follows from rig (, 2 d < 1 where 

n^ d ^ is the crossing point of 774 and 775. This yields the 
requirement 



£ < a 1 / 6 ^) 1 / 2 , (d<d*). 



(45) 



The crossover from widely to narrowly spaced layers can 
be obtained from the criterion Z s = d. It occurs at the 
spacing 



7 j* "-z -2/3 

d = d = — a ~ a ' a. 

K 



In the widely spaced case one hnds 



e 2 N' 



,V2 



A< 



K 3 / 2 rid> K 3 / 2 2 71(2' 

and thus the corrected transition line 



714 



/*,y/4 AT 1 / 2 _ _ 1/e N^_ 



(46) 



(47) 



(48) 



The Mott transition at low density now takes place at 
a density 



-J 2d ) - 



riM = n B 



1 _ a 1 ' 3 
i 2 <i a 2 d ' 



(49) 



where a x = a _1 / 6 a is the effective Bohr radius in 
the plane. It is obtained by comparing kinetic and 
Coulomb energy in the plane, ft 2 / mtt i = ^ l\^^ a 'x\ = 
e 2 / '[k 3 / 2 (a x I Kx )]■ Logarithmic corrections occur in this 
case if 77m<2£ 2 < 1, i.e., for 



( < oq- 1 / 6 , (d>d*). 



(50) 



Note that the critical lines 714 for narrowly and widely 
spaced systems match when d — d* . The same holds 
for the critical size of pairs necessary to have a BEC-like 
regime. 



C. Is the BEC limit of this theory applicable to 
high T c superconductors? 

Above we have obtained results along two lines. First, 
for relatively large pairs (in the BCS regime) we predict 
the SI transition line given by Eqs. (|3"6")) and (|57|) . Sec- 
ond, for small pairs we have discussed an additional log- 
arithmic factor originating from BEC effects. One may 
question whether the value of £ in high T c superconduc- 
tors is small enough so that the condition nd£ 2 < 1 for a 



BEC-like regime with extra logarithmic factors is realis- 
tic. The most frequently cited number for under-doped 
uncompensated YBCO for the superconducting coher- 
ence length is 2 nm. However, for the case of strong 
coupling of holes the size of pairs, £, can be smaller than 
the coherence length. ARPES data in YBCO indicate 
that actually 33 £ ~ lnm. Using that that the bound- 
ary of the superconducting dome on the underdoped side 
occurs at (n u d) — 0.06<2q 2 where a,Q ~ 0.4nm is the lat- 
tice constant of the two-dimensional Cu lattice, one hnds 
that the SI transition at T = empirically happens when 
nd£ 2 ~ 0.4 < 1. This means that the BEC part of our 
diagram Fig. [4] is marginally relevant. In iron- arsenide 
superconductors the small value £ < 2 nm was recently 
found 3 -^. Since n u da^ has a similar value as in cuprates, 
this leads to n u d^ 2 ~ 1 in this new family of supercon- 
ductors as well. 

The low density (BEC) regime might indeed be exper- 
imentally relevant if one adopts a popular interpretation 
of the pseudogap which is observed in underdoped sam- 
pj es 23 i 35 l 36 i 37 i 38 i 39^ j a tt e r assumes that the pseudo- 
gap is due to preformed hole pairs with large binding 
energy (-E pa ir > T c ), the pairs being localized by disor- 
der at low doping density. If such an interpretation is 
correct, the small A^ part of our diagram Fig. 0]may be 
relevant for high T c superconductors. 

As we have seen the SI boundary reflects the BEC- 
BCS crossover in the form of extra logarithmic factors in 
the critical density 714 only if £ is sufficiently small, i.e., 
if condition or (|50"|) is satisfied. Let us discuss this 
condition for the example of Bi2Sr2Cu206+5 (Bi-2201). 
The mean distance between copper planes is d — 12. 3A, 
and the lattice spacing in the planes is ao — 5.36A. From 
recent optical measurements^, one can extract the effec- 
tive mass of carriers as m e a ~ 3 — 4m e in the under- 
doped regime. The dielectric constant along the c-axis 
isi 2 . k z = 18.9. We are not aware of direct measurements 
of k x , but usually, the anisotropy is relatively modest^, 
e.g., a = k x /k z = 0.99 in Nd 2 Cu0 4 , or 0.7 in Pr 2 Cu0 4 . 
Neglecting the anisotropy we can use k — k z and the 
effective mass to estimate the Bohr radius in Bi-2201 as 
a ~ a x ~ a z rs 3 A. We can compare this to an al- 
ternative estimate obtained as follows: We assume that 
the SI transition of uncompensated materials is essen- 
tially a Mott transition of doped carriers, which is known 
to occur roughly wheni 3 - (n M d)a 2 = c with^ c w 0.04. 
Using an approximate value for n u d « 0.06etQ 2 we find 
a w 4.4 A, in rough agreement with the above calculation 
based on the effective mass. This material thus certainly 
corresponds to widely spaced layers, d > a. The esti- 
mated Bohr radius a is of the same order as the typical 
pair size £ in strongly underdoped samples. The require- 
ment ([50")) for observing the SI transition in the BEC 
regime is thus just marginally satisfied in this standard 
cuprate compound. 

More favorable conditions for the crossover to the BEC 
limit may be expected in materials with high dielectric 
constants (such as in La2- a; Sr a ;Cu04), which increases 
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the Bohr radius. A similar tendency can be expected 
from a small effective mass (small band mass and/or 
small mass renormalization) , provided it does not occur 
simultaneously with an increase of the pair size £. 

The application of our theory to certain particular su- 
perconductors may require further adjustments of the 
model. For example, in Yi_ z La z (Bai_ x La x )2Cu30 y ac- 
ceptors are divalent and we have to define proper vari- 
ables for the phase diagram. We can use Na — (y— 6) /v uc 
for the concentration of divalent oxygen acceptors (exces- 
sive oxygen) , Njy — x /v uc for the concentration of mono- 
valent donors, and n — 2N A ~ Np for the concentration 
of holes. Here v uc is the volume of the unit cell. It is easy 
to show that the concentration N — 4N A + No plays the 
role of the effective concentration of monovalent charged 
impurities. Indeed, for randomly distributed impurities 
in a given volume R 3 , the variances of the donor, ac- 
ceptor and net charge number distribution are equal to 
N A R 3 , N D R 3 and (p 2 N A + N D )R 3 respectively, where 
p represents the valence of the acceptor. So the effec- 
tive concentration of monovalent charged impurities is 
not N — N A + but N = p 2 N A + N D . For excess 
oxygen atoms one has p = 2, and the coefficient 4 in the 
expression for N reflects the enhanced role of divalent 
charge in the creation of potential fluctuations. 

There is a further complication for YBCO, in that a 
fraction of holes does not reside in CuO planes, but in 
CuO chains. This should be taken into account when 
comparing our theory with YBCO data. However, most 
other high T c layered superconductors do not suffer from 
such a complication. 



IV. CONCLUSION 

In conclusion we have established the phase diagram 
for the superconductor-insulator transition in heavily 
doped, strongly compensated semiconductors endowed 
with a strong superconductive coupling mechanism. The 
phase transition line at large impurity and carrier density 
coincides essentially with the well-known metal-insulator 
transition in doped semiconductors. However, if Cooper 
pairs are tightly bound, such that £ < a, there is a low 
density (BEC) regime where preformed pairs are dilute 
even at the SI transition. In this regime we have estab- 
lished two new segments of the SI transition line which 
reflect that a gas of compact bosons is more compressible 
than an equally dense gas of weakly interacting fermions. 

Recently, an interesting system exhibiting a direct SI 
transition upon doping has been discovered in the form 
of boron-doped diamond^. The latter can be simulta- 
neously doped by both donors and acceptors and thus 
constitutes a promising system in which one might ob- 
serve the effects we predict for the isotropic 3d case. It 
would also be interesting to test our predictions numeri- 
cally, e.g., following the lines of recent work which investi- 
gated the interplay of superconductivity and localization 
in strong disorder—. 



We have extended considerations from the isotropic 
case to layered systems such as the cuprates. In this case, 
we have found new equations for the SI line which can 
be verified experimentally. We showed that due to the 
smaller phase space for bosons in the plane, the crossover 
to the BEC regime manifests itself on the phase transi- 
tion line only by an additional logarithmic factor. Ex- 
cept for the logarithmic factors, our results also apply to 
the metal- insulator transition in layered, strongly doped 
fermion systems. 

Apart from determining the phase transition line, we 
have established the properties of the insulating phase. 
We assert that in the presence of strong superconduct- 
ing couplings all fermions are paired, and hence at low- 
est temperatures, transport is due to the variable range 
hopping of Cooper pairs. This kind of transport may be 
observable in compensated diamond and high T c materi- 
als. 
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APPENDIX A: SMALL PAIRS IN BIPOLARON 
SYSTEMS 

Frohlich polarons are quasiparticles arising in systems 
with strong electron phonon coupling. An extensive 
study of polarons is given in Ref. [47], based on Feynman's 
path integral approach, giving accurate results in dimen- 
sions d = 2,3. If the coupling strength is large enough 
polarons can bind into strongly bound pairs which finally 
undergo a SI transition if they are dense enough. 

An essential ingredient for strong coupling is a signif- 
icantly small ratio between the electronic dielectric con- 
stant K e i, and the static one, n > K e i, 



T)= . 

K 



(Al) 



Note that the static dielectric constant is the one which 
enters the screening problems discussed in the main text. 

The electron-phonon coupling is characterized by the 
coupling constant, 



e 2 / 2mu LO \ 1/2 
2htOLO \ h, J 



1 

Kel 



1 



(A2) 



e 2 / 2mujLo\ 
2fiuj L o \ K J 



1/2 



— = — ■ (A3) 

^el ^oo 



Here, m is the band mass and u>lo the long wavelength 
optical phonon frequency. In the last step it was assumed 
that i] <C 1. Further, 



me 2 



< a, 



(A4) 
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is an "effective Bohr radius" built from the electronic 
dielectric constant. However, the Bohr radius which ap- 
pears in the theory presented in the main text is given 
by 



h 2 K 

n po \e- 



(A5) 



where m po i > m is the polaron mass. At large coupling, 
one finds 4 ^ m po i/m oc a. Note that a can be much larger 
than Ooo if r\a <C 1. 

At strong coupling polarons can bind into pairs, so- 
called bipolarons. The bipolaron radius is usually only 
10 — 20% larger 47 than the radius of a single polaron, 
which is given by 



muj L o 



1/2 



(A6) 



These pairs are stable if the coupling is sufficiently 
strong (a > a c ), and if the ratio rj is sufficiently small. 
The critical values have been computed 47 to be a c = 2.9 
(2D) and 6.8 (3D). 

For not too strong couplings a > a c the maximal 
admissible ratio of dielectric constants which allows for 
bipolaron formation was found to be (for both 2d and 
3d) 



(7 - l)(q - q c ) 
j(a - a c ) + a c 



(A7) 



with 7 w 1.63/V2 w 1.15. Note that r/ c -> as a — » a c . 

We are eventually interested in the possibility of small 
pairs with £ < a. Note that even though a large coupling 
a implies £ > , a small ratio 77 <gc 1 still allows one to 
have 



£ m po\ 2 . , 

— = a 77 — — ~ a 77 < 1 . 

a m 



(A8) 



This is precisely what is needed to observe the BEC part 
of the SI transition line. 

At very strong coupling, larger conglomerates of po- 
larons can form stable bound states. The formation of 
such multipolarons has not been studied systematically 
yet, apart from an analysis at asymptotically strong cou- 
pling 49 . For multipolarons to exist, the coupling needs 
to be considerably stronger, a 3> a c . At a given large 
a, the most stable bound state will depend on the ratio 
of dielectric constants, 77. The smaller 77, the more po- 
larons can bind together. For example, at asymptotically 
strong coupling one finds bipolarons for 0.046 = rj c < 
t] < T]c — 0.079, and larger multipolarons for 77 < r/^ . 

Multipolarons containing an even number of polarons 
will be compact bosons, which eventually undergo a Bose 
Einstein condensation in sufficiently weak disorder. Up 
to numerical prefactors, the Si-transition for such multi- 
polarons would be of the same nature as the one discussed 
in the main text for bosons formed by pairs of carriers. 
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